Riemann–Siegel formula

In mathematics, the Riemann–Siegel formula is an asymptotic formula for the error of the approximate functional equation of the Riemann zeta function, an approximation of the zeta function by a sum of two finite Dirichlet series. It was found by Siegel (1932) in unpublished manuscripts of Bernhard Riemann dating from the 1850s. Siegel derived it from the Riemann–Siegel integral formula, an expression for the zeta function involving contour integrals. It is often used to compute values of the Riemann–Siegel formula, sometimes in combination with the Odlyzko–Schönhage algorithm which speeds it up considerably. When used along the critical line, it is often useful to use it in a form where it becomes a formula for the Z function.

If M and N are non-negative integers, then the zeta function is equal to

$\zeta(s) = \sum_{n=1}^N\frac{1}{n^s} %2B \gamma(1-s)\sum_{n=1}^M\frac{1}{n^{1-s}} %2B R(s)$

where

$\displaystyle\gamma(s) = \pi^{1/2-s}\Gamma(s/2)/\Gamma((1-s)/2)$

is the factor appearing in the functional equation ζ(s) = γ(s) ζ(1 − s), and

$R(s) = \frac{-\Gamma(1-s)}{2\pi i}\int \frac{(-x)^{s-1}e^{-Nx}dx}{e^x-1}$

is a contour integral whose contour starts and ends at +∞ and circles the singularities of absolute value at most 2πM. The approximate functional equation gives an estimate for the size of the error term. Siegel (1932) and Edwards (1974) derive the Riemann–Siegel formula from this by applying the method of steepest descent to this integral to give an asymptotic expansion for the error term R(s) as a series of negative powers of Im(s). In applications s is usually on the critical line, and the positive integers M and N are chosen to be about (2π Im(s))^1/2. Gabcke (1979) found good bounds for the error of the Riemann–Siegel formula.

Riemann's integral formula

Riemann showed that

$\int_{0\searrow 1}\frac{e^{-i\pi u^2%2B2\pi i pu}du}{e^{\pi i u}-e^{-\pi i u}} = \frac{e^{i\pi p^2}-e^{i\pi p}}{e^{i\pi p}- e^{-i\pi p}}$

where the contour of integration is a line of slope −1 passing between 0 and 1 (Edwards 1974, 7.9).

He used this to give the following integral formula for the zeta function:

$\displaystyle\pi^{-s/2}\Gamma(s/2)\zeta(s)=$

$\pi^{-s/2}\Gamma(s/2)\int_{0\swarrow 1}\frac{x^{-s}e^{\pi i x^2}}{e^{\pi i x}-e^{-\pi i x}}\,dx %2B\pi^{-(1-s)/2}\Gamma((1-s)/2)\int_{0\searrow 1}\frac{x^{s-1}e^{-\pi i x^2}}{e^{\pi i x}-e^{-\pi i x}}\,dx$

References

Berry, Michael V. (1995), "The Riemann–Siegel expansion for the zeta function: high orders and remainders", Proceedings of the Royal Society. London. Series A. Mathematical, Physical and Engineering Sciences 450 (1939): 439–462, doi:10.1098/rspa.1995.0093, ISSN 0962-8444, MR 1349513

Edwards, H. M. (1974), Riemann's Zeta Function, New York: Dover Publications, ISBN 978-0-486-41740-0, MR 0466039
Gabcke, Wolfgang (1979), Neue Herleitung und Explizite Restabschätzung der Riemann-Siegel-Formel, Georg-August-Universität Göttingen, http://www.worldcat.org/oclc/60563505
Siegel, C. L. (1932), "Über Riemanns Nachlaß zur analytischen Zahlentheorie", Quellen Studien zur Geschichte der Math. Astron. und Phys. Abt. B: Studien 2: 45–80 Reprinted in Gesammelte Abhandlungen, Vol. 1. Berlin: Springer-Verlag, 1966.

Riemann–Siegel formula

Riemann's integral formula

References

External links